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ABSTRACT 

We calculate the photonuclear production of heavy quarks in ultraperipheral heavy ion 
collisions. The integrated cross section and the rapidity distribution are computed employing 
sound high energy QCD formalisms as the collinear and semihard approaches as well as the 
saturation model. In particular, the color glass condensate (CGC) formalism is also considered 
^ . using a simple phenomenological parameterization for the color field correlator in the medium, 
which allow us to obtain more reliable estimates for charm and bottom production at LHC 
^ \ energies. 
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o : 1 Introduction 

The recent results from RHIC suggest that relativistic heavy ion collisions at high energies 
probe QCD in the non-linear regime of high parton density [1]. In such a regime the growth 
of parton distributions should saturate, possibly forming a color glass condensate [2] (for a 
pedagogical presentation see Refs. [3, 4]), which is characterized by a bulk momentum scale 
^ ■ Q s . In particular, the RHIC data on the multiplicity distribution of the produced hadrons 
as a function of centrality, rapidity, and collision energy appears to be consistent with the 
CGC predictions. However, there are still a number of open questions, mainly associated to 
the fact that other models based on different assumptions reasonably describe the same set of 
experimental data [5, 6]. For instance, the recent RHIC data indicate the lack of suppression in 
the high-pr spectra of charged hadrons produced in d + Au collisions [7], in contrast with the 
initial expectation of the CGC formalism [1]. The main uncertainty present in those analysis is 
directly connected with the poor knowledge of the initial conditions of the heavy ion collisions. 
Theoretically, the early evolution of these nuclear collisions is governed by the dominant role 
of gluons [8], due to their large interaction probability and the large gluonic component in the 
initial nuclear wave functions. Consequently, a systematic measurement of the nuclear gluon 
distribution is of fundamental interest to understand the parton structure of nuclei, determine 
the initial conditions of the QGP and constraint the QCD dynamics at high energies. 



In the last years, there has been a lot of interest in the description of electron-nucleus 
collisions at high energies, with particular emphasis in the behavior of the nuclear structure 
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functions [9] and their logarithmic slope [10] at small value of the Bjorken x variable, obtaining 
predictions which agree with the scarce experimental data. Moreover, the high energy heavy 
quark photoproduction on nuclei targets has been studied in detail in Refs. [11, 12], considering 
the available several scenarios for the QCD dynamics at high energies. The results of those 
analysis show that future electron-nucleus colliders at HERA and RHIC [13, 14], probably could 
determine whether parton distributions saturate and constraint the behavior of the nuclear 
gluon distribution in the full kinematical range. However, until these colliders become reality 
we need to consider alternative searches in the current and/or scheduled accelerators which allow 
us to constraint the QCD dynamics. Here, we analyze the possibility of using ultraperipheral 
heavy ion collisions as a photonuclear collider and study the heavy quark production assuming 
distinct approaches for the QCD evolution. 

The analysis of heavy quark production in ultraperipheral heavy ion collisions has been 
proposed many years ago [15], improved in Refs. [16, 17, 18] and recently revisited in Refs. 
[19, 20, 21]. The most of these approaches calculates the cross section assuming the validity 
of the collinear factorization, where the cross sections involving incoming hadrons are given, 
at all orders, by the convolution of intrinsically non-perturbative, but universal, quantities - 
the parton densities, with perturbatively calculable hard matrix elements, which are process 
dependent. In this approach all partons involved are assumed to be on mass shell, carrying only 
longitudinal momenta, and their transverse momenta are neglected in the QCD matrix elements. 
However, in the large energy (small-x) limit, the effects of the finite transverse momenta of 
the incoming partons become important, and the factorization must be generalized, implying 
that the cross sections are now /cj_-factorized into an off-shell partonic cross section and a k±- 
unintegrated parton density function J r (x,k±), characterizing the /cj_-factorization approach 
[22, 23, 24]. The function T is obtained as a solution of the evolution equations associated to 
the dynamics that governs the QCD at high energies. Here, we estimate, for the first time, 
the total cross section and the rapidity dependence of the nuclear photoproduction of heavy 
quarks in ultraperipheral heavy ion collisions, considering the /c^-factorization approach and 
distinct nuclear unintegrated gluon distributions. Moreover, we extend the previous study for 
heavy quark production in the color glass condensate formalism considering a realistic photon 
flux and a phenomeno logical dipole cross section which is energy dependent, allowing to obtain 
reliable estimates for the rapidity distribution within this formalism. For comparison, we also 
present the predictions for the cross section from the collinear factorization approach. 

This paper is organized as follows. In next section we present a brief review of ultraperipheral 
heavy ion collisions and write down the main formulas describing the photonuclear process in 
these collisions. In Section 3 we discuss some sound models for the heavy quark photoproduction 
in the collinear and /^-factorization approaches. Moreover, we discuss our main assumptions 
in order to extend the previous results on color glass condensate formalism. Finally, in Section 
4 we present our results for the total cross section and rapididity distribution of charm and 
bottom production for the LHC energies. 
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2 Ultraperipheral relativistic heavy ion collisions 



In heavy ion collisions the large number of photons coming from one of the colliding nuclei will 
allow to study photoproduction, with energies W^n reaching to 950 GeV for the LHC. The 
photonuclear cross sections are given by the convolution between the photon flux from one of 
the nuclei and the cross section for the scattering photon-nuclei. The photon flux is given by 
the Weizsacker- Williams method [25] . The flux from a charge Z nucleus a distance b away is 



d 3 N(uj, b 2 ) Z 2 a ernV 
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where 7^ is the Lorentz boost of a single beam and rj = oub/^L', Koiv) an d K\{rj) are the mod- 
ified Bessel functions. The requirement that photoproduction is not accompanied by hadronic 
interaction (ultraperipheral collision) can be done by restricting the impact parameter b to be 
larger than twice the nuclear radius, Ra = 1.2 A 1 / 3 fm. Therefore, the total photon flux inter- 
acting with the target nucleus is given by Eq. (1) integrated over the transverse area of the 
target for all impact parameters subject to the constraint that the two nuclei do not interact 
hadronically. An analytic approximation for AA collisions can be obtained using as integration 
limit b > 2 Ra, producing 
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where 77 = 2uj Ra/^l- The final expression for the production of heavy quarks in ultraperipheral 
heavy ion collisions is then given by, 

00 dN ( ) 

a AA^QQX (v^Nn) = / duJ duJ a 1 A^QQx{ W ^A = 2uJ \fSN^j ( 3 ) 

where uj m in = ^Qq/^lL'mp and V^nn is the c.m.s energy of the nucleus-nucleus system. The 
Lorentz factor for LHC is 7^ = 2930, giving the maximum c.m.s. 7iV energy W 1 a ~ 950 GeV. 
It is worth mentioning that the difference between the complete numeric and the analytical 
calculation presented above for the photon flux is less than 15 % for the most of the purposes 
[25]. 

Before considering the distinct models for the photon-nucleus cross section, it is interesting 
to determine the values of x which will be probed in ultraperipheral heavy ion collisions. The 
Bjorken x variable is given by x = (Mgg/2p)e~^, where Mqq is the invariant mass of the 
photon-gluon system and y the center of momentum rapidity. For Pb + Pb collisions at LHC 
energies the nucleon momentum is equal to p — 2750 GeV; hence x = (Mqq/ 5500 GeV) e~ y . 
Therefore, the region of small mass and large rapidities probes directly the high energy (small 
x) behavior of the QCD dynamics present in the 7 A cross section. This demonstrates that 
ultraperipheral heavy ion collisions at LHC represents a very good tool to constraint the high 
energy regime of the QCD dynamics, as already verified for two-photon processes [26, 27, 28]. 
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3 Models for nuclear heavy quark production 



For our further analysis on photonuclear production of heavy quarks we will consider distinct 
available high energy approaches. First, we take into account the usual collinear approach, 
where the production cross section is driven by the collinear gluon distribution on the nuclei. 
This one contains a lot of information about nuclear shadowing, EMC, and anti-shadowing 
effects. Second, the fc^-factorization formalism is introduced, where the relevant quantity is 
now the nuclear unintegrated gluon distribution. For this purpose, we analyze two simple 
parameterization for it which are consistent with data description on inclusive and diffractive 
DIS. Finally, we take into account the color glass condensate formalism, where the scattering 
process is viewed as the interaction of the probe particles with the strong nuclear color field 
treated in a classical approximation. The main quantity is the correlator of two Wilson lines, 
which is related to the dipole cross section and at lowest order has no energy dependence. We 
present a simple phenomenological parameterization which introduces higher orders corrections 
to the classical approximation and allow us produce more realistic estimates for the cross section. 

3.1 The collinear approach 

In hard photon-hadron interactions the photon can behave as a pointlike particle in the so- 
called direct photon processes or it can act as a source of partons, which then scatter against 
partons in the hadron, in the resolved photon processes (For a recent review see Ref. [29]). 
Resolved interactions stem from the photon fluctuation into a quark-antiquark state or a more 
complex partonic state, which are embedded in the definition of the photon structure functions. 
Recently, the contribution of resolved photon processes in ultraperipheral heavy ion collisions 
was discussed in Ref. [20] and studied in detail in Ref. [21]. One of the main results is that at 
LHC, these contributions are ~ 15 and 20 % of the total charm and bottom photoproduction 
cross sections, respectively, comparable to the shadowing effect. Here, we will consider only the 
direct photon contribution. 

At high energies the main subprocess occurring when the photon probes the structure of 
the nucleus is the photon-gluon fusion producing the heavy quark pair. It can be described 
through perturbative QCD, with the cross section given in terms of the convolution between 
the elementary cross section for the subprocess — > QQ and the probability of finding a gluon 
inside the nucleus, namely the nuclear gluon distribution. In this collinear approach the heavy 
quark photoproduction cross section is given by 

fW 2 , da — 

^ Q q(W,a) = J I dM^-p^g A (x,S), (4) 

J4m Q U1V1 QQ 

where the quantity dagg/dM^ is calculable perturbatively, Mqq is the invariant mass of the 
heavy quark pair with x = M^/W 2 A . The c.m.s energy of the system is labelled W lA and 

g A (x,fi 2 ) is the gluon density inside the nuclear medium, with \x being the factorization scale 
and rriQ the heavy quark mass. For our purpose here we will use fj 2 = 4mg, with m c = 1.5 
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GeV and m b = 4.5 GeV. The differential cross section in leading order is given by [30] 



M 2 - 

QQ 



(5) 



where 6q is the heavy quark charge and (3 = Arriq/M 2 -^. In our further calculation on the 

collinear approach one considers that xqa{x,Q 2 ) = R g (x,Q 2 ) x xgN{x,Q 2 ), where R g param- 
eterize the medium effects as proposed in Ref. [31] and xgj^ is the nucleon gluon distribution 
given by the GRV98(LO) parameterization [32]. It is worth mentioning that different choices 
for the factorization scale and quark mass produce distinct overall normalization to the total 
cross section at photon-nucleus and ultraperipheral nucleus-nucleus interactions. For details 
see Ref. [12] , where the heavy quark photoproduction at eRHIC and THERA energies has been 
discussed. In the Section 4 we discuss the dependence of our results on the choice of the quark 
mass and parton distribution parameterization. 



3.2 The fc^-factorization formalism 

In the /c^-factorization (or semihard) approach, the relevant QCD diagrams are considered with 
the virtualities and polarizations of the initial partons, carrying information on their transverse 
momenta. The scattering processes are described through the convolution of off-shell matrix 
elements with the unintegrated parton distribution, J-(x, k±) (see [33] for a review). The latter 
can recover the usual parton distributions in the double logarithmic limit by its integration 
over the transverse momentum of the k± exchanged gluon. The gluon longitudinal momentum 
fraction is related to the c.m.s. energy, W 7 a, in the heavy quark photoproduction case as 
x = 4m 2 Q/W 2 A) as in the collinear case. The cross section for the heavy-quark photoproduction 
process is given by the convolution of the unintegrated gluon function with the off-shell matrix 
elements [33, 34, 35, 36]. Considering only the direct component of the photon we have that 
°foT* reads as [35], 

phot/rrr \ ft ™ e Q f j ,2 j2i a s{l^ 2 ) ^(x, k\; ii 2 ) 
(hot \Wia) = / dz d p 1± d k± — 2 

where D\ = p 2 ± + rrig and D 2 = (k± — p 1A _) 2 + nig. The transverse momenta of the heavy 
quark (antiquark) are denoted by p 1± and p 21 _ = (k± — p u ), respectively. The heavy quark 
longitudinal momentum fraction is labeled by z. The scale fi in the strong coupling constant in 
general is taken to be equal to the gluon virtuality, in close connection with the BLM scheme 
[37]. Here, we will use the prescription /i 2 = k 2 ± + m 2 Q. 

In order to perform a phenomenological analysis within the /c^-factorization approach, in 
the following we use two distinct parameterizations for the unintegrated gluon distribution (For 
details see Ref. [12]). First, one considers the derivative of the collinear nuclear gluon parton 
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distribution function, quite successful in the proton case and tested in the nuclear case in Ref. 
[12]. It simply reads as, 



T ( r 1,2 . A \ _ dxG A (x, k\) 
•/ nuc \ x , K ±) ^) — d\nk 2 1 

where xGa(x, Q 2 ) is the nuclear gluon distribution, which was taken from the EKS parame- 
terization [31] for the medium effects and the GRV94(LO) for the nucleon parton distribution 
[38]. The latter choice is supported by the good description of heavy quark photoproduction 
in the full kinematical region [12]. As a consequence, with this procedure we include in our 
calculations the medium effects (shadowing, antishadowing, EMC and Fermi motion effects) 
estimated by that parameterization. Moreover, we emphasize that this nuclear gluon distribu- 
tion is solution of the DGLAP evolution equations, which is associated to a linear dynamics 
that does not consider dynamical saturation effects. 

The second parameterization is given by the model introduced in Ref. [39], which provides 
an extension of the ep saturation model through Glauber-Gribov formalism. In this model the 
cross section for the heavy quark photoproduction on nuclei targets is given by [39, 11] 

alf(W,A) =J Q 1 dz J d 2 r \* T (z, r, Q 2 = 0)| 2 aj ip (x, r 2 ,A) , (8) 

where the transverse wave function is known (See e.g. Ref. [40]). As \^l\ 2 oc Q 2 , the longitu- 
dinal piece does not contribute for Q 2 = 0. The nuclear dipole cross section is given by, 

r *> A ) = j d% 2 { 1 - ex P [-\ AT A(b) a^{x, r 2 )] } , (9) 

where b is the impact parameter of the center of the dipole relative to the center of the nucleus 
and the integrand gives the total dipole-nucleus cross section for fixed impact parameter. The 
nuclear profile function is labelled by T^b), which will be obtained from the 3-parameter 
Fermi distribution for the nuclear density [41]. The parameterization for the dipole cross 
section takes the eikonal-like form, a^ ip (x, r 2 ) = <7o[l — exp (— Q 2 (x) r 2 /4) ], where one has 
used the parameters from [42], which include the charm quark with mass m c = 1.5 GeV and 
the definition x = (Q 2 + Am 2 Q)/W 2 A . The saturation scale Q 2 (x) = (x /x) x GeV 2 , gives the 
onset of the saturation phenomenon to the process. 

The equation above sums up all the multiple elastic rescattering diagrams of the qq pair 
and is justified for large coherence length, where the transverse separation r of partons in the 
multiparton Fock state of the photon becomes as good a conserved quantity as the angular 
momentum, namely the size of the pair r becomes eigenvalue of the scattering matrix. The 
corresponding unintegrated gluon distribution can be recovered from a Bessel- Fourier transform 
to the momentum representation 



= B^fH^ izgexp (10) 

7T 2 a s \QtJ m =ln=0 ml 11 \ TlQjJ 
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which depends on the transverse momentum k± through the scaling variable r = k\/Q 2 . The 
unintegrated gluon vanishes asymptotically at k\ — > 0, oo and its maximum can be identified 
with the saturation scale Q s a(x) [39, 43]. The model recovers the original one for the proton 
case, taking A = 1 and the normalization condition / d 2 bTA(b) = 1. 

In Ref. [11], it has verified that the resummation at the proton level is less sizeable in 
the final results at nuclear level. Therefore, this fact allow us take just the color transparency 
behavior on the dipole nucleon cross section. Hence, in such a particular case one can compute 
analytically the unintegrated gluon distribution, which is expressed as 

= ^ (^_) -p(-q^)) • (11) 

where Q 2 A (x) = ^AT A (b) <JqQ 2 {x) define the nuclear saturation scale. Such an approximation 
is justified in the heavy quark case, which is dominated by small dipole configurations (large 
transverse momentum k\ ~ mq). It is clear that Eq. (11) presents a scaling pattern on 
the variable r = k 2 ± /Q 2 A , which implies scaling on r in the nuclear heavy quark production. 
Recently, this feature has been shown also in the nucleon case [44]. 

Some comments are in order here. In our analysis we have disregarded higher-order Fock 
states in the photon wave function and considered only the evolution in the dipole cross section. 
In principle, at leading order and for inclusive process this is a reasonable approximation. 
However, for diffractive processes, for example, the qqg component cannot be disregarded. The 
recent results for the NLO corrections for the impact factor [45] will allow to verify the validity 
of the color dipole approach at higher orders. Other point which deserves discussion is that 
we have assumed the validity of the /c_|_-factorization in photon-nucleus interactions. At high 
energies and intermediate densities, it is a reasonable assumption, since the derivation of the 
^-factorization presented in Ref. [46] can be directly extended for the nuclear heavy quark 
photoproduction. However, for very large parton densities a breakdown of the k± (and collinear) 
factorization is expected, mainly associated to the effects of the non-linearity of the non-Abelian 
gluon field [2, 47, 48]. Finally, we also have disregarded the resolved photon contribution in the 
semihard approach. For completeness, lets quote the estimates for the resolved component on 
the nucleon level. In Ref. [36] it gives a contribution of order 20-30 %, rising faster on energy 
(hence important mostly at large rapidities) than the direct photon component and being stable 
under different choices for the unintegrated gluon distributions and quark mass. Moreover, the 
resolved contribution is somewhat higher than in the collinear approach due to the non-zero 
transverse momentum transfer effect. We expect a similar trend in photon-nucleus interactions. 



3.3 The Color Glass Condensate formalism 

At small x and/or large A one expects the transition of the regime described by the linear 
dynamics (DGLAP, BFKL) (For a review, see e.g. Ref. [49]), where only the parton emissions 
are considered, to a new regime where the physical process of recombination of partons becomes 
important in the parton cascade and the evolution is given by a nonlinear evolution equation. In 
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this regime a Color Glass Condensate (CGC) is expected to be formed [2], being characterized 
by the limitation on the maximum phase-space parton density that can be reached in the 
hadron/nuclear wavefunction (parton saturation) and very high values of the QCD field strength 
F[iv ~ [50]. The large values of the gluon distribution at saturation (large occupation 

number of the soft gluon modes) suggests the use of semi-classical methods, which allow to 
describe the small- x gluons inside a fast moving nucleus by a classical color field. This color 
field is driven by a classical Yang-Mills equation whose source term is provided by faster partons. 
When the energy further increases, the structure of the classical field equations does not change, 
but only the correlations of the color source. This change can be computed in perturbation 
theory, and expressed as a functional renormalization group equation for the weight function, 
in which the 'fast' partons are integrated out in steps of rapidity and in the background of 
the classical field generated at the previous step. This approach enables one to calculate cross 
section cross sections in a high gluon density environment. Recently, this approach has been 
applied for eA [51, 47], pA [47, 52, 53] and AA processes [48]. 

In Refs. [19] the heavy quark production in ultraperipheral heavy ion collisions has been 
analyzed in the color glass condensate formalism. In particular, those authors have considered 
the photon-nuclei interaction, taking into account the electromagnetic interaction to lowest 
order in the coupling and the interactions with the strong color background field to all orders. 
The quantum evolution is not included in the calculations. Their prediction for the rapidity 
(y) distribution of the heavy quark (or antiquark) is given by 



TL 



<** - nR l j^gOlj 7 ft 7 ik l 0( » ±) [ 1 + 4{kl - ZD arcth 



dy for* ^ b< J >\ k ± ^T^ jul + 4m 2 Q 

with N c being the color number and e q the quark charge. The color field correlator C(k±) in 
the medium is given by 

C(k ± ) = Jd 2 re tk ^ r e- B ^ = Jd 2 re lk ^ r (u{0)tf(r)) (13) 

with (..) representing the average over all the configurations of the color fields in the nucleus. 
The unitarity matrix U(r) contains the information related to the interactions between the 
quark and the colored glass condensate (classical color field of the nucleus) and is expressed in 
terms of the color sources in the nucleus. Therefore, C(k±) depends on the structure of the 
color sources describing the target nucleus and describes the interactions of a high energy probe 
with the target. 

In Ref. [19] the McLerran-Venugopalan model for the correlator C(k±) was considered. In 
this model the function B 2 (r) in Eq. (13) is approximated as follows 

B 2 (r)*9kll ln (_L) , (14) 



47r \r A 



where the saturation scale Q s a at this classical level does not depend on the rapidity (energy) 
and A is an infrared cutoff related to the scale at which color neutrality occurs. In principle, 



A is at least as large as Aqqd- For a saturated target, it can be probed that color neutrality 
occurs over transverse spatial scales as small as 1/Q S [54]. Lets summarize the results coming 
from applying the approximation above and A = Aq C d m Eq. (12). Supposing Q s a/Aqcd = 10 
(y = 2.3) one obtains da/dy = 355 mb for charm and da/dy = 11 mb for bottom with the 
rapidity distribution being flat on y. 

In principle, extending the previous calculation to include quantum evolution is just a mat- 
ter of changing the C(k±) present in the previous calculations by one that has been calculated 
considering the evolution of the sources. However, the general solution of the functional renor- 
malization group equation for the weight function is not known, but only approximate solutions 
in some limiting kinematical regimes have been obtained [55]. In general, the quantum correc- 
tions lead to a modification of the distribution function of the hard sources when the energy 
increases, which implies a rapidity dependence for the correlator C(k±). 

In order to go further and introduce an energy dependence on the calculations one makes 
use of the fact that the function C(k±) is directly related to the Fourier transform of the 
dipole-nucleus total cross section, as follows 

C (x, k_i) = -^-^ / d2r e 1 k± r Wdi P (x,r -> oo) -a dip (x,r)} . (15) 

For this definition of the correlator, we have a direct relation between this quantity and the 
unintegrated gluon density, given by J-(x, k±) = (3R\/8ir 2 a s ) k 2 ± C (x, k±). In our analysis we 
will employ an educated guess for the dipole-nucleus cross section 



a£ p (x,r) = 27rR 2 A 



1 — exp 



(16) 



where we assume Q 2 A (x) = A 1 / 3 x Q 2 (x). This model is inspired in the saturation model for 
the dipole-proton scattering proposed in Ref. [42], and encodes the main properties of the 
high density approaches, namely color transparency for small pair separations (r — > 0) and 
saturation for large pair separations (r — > oo). Moreover, the nuclear A-dependence assumed 
for the saturation scale agrees with the recent analysis for the scattering of a small dipole in 
a nuclear target in the McLerran-Venugopalan model and fixed coupling BFKL dynamics [56]. 
However, the energy dependence of the saturation scale still is an open question. Recently, 
several groups have studied the solutions of the Balitski-Kovchegov equation [57] (See also Ref. 
[58]), with x dependence for Q s slightly different from the ansatz used here. 

Embeding Eq. (16) in the definition on Eq. (15) we obtain the following analytical expres- 
sion for the color field correlator 



e(l,fe±) = e^ exp (~efe))' (17) 



which also obeys the normalization condition / d 2 kC(x, k±)/(2ir) 2 = I. In the next section, we 
compute the photonuclear cross section using the correlator above and in addition consider the 
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Figure V. The rabidity distribution for the distinct high energy approaches: collinear approach 
(long-dashed lines), semihard formalism (solid and dotted lines), saturation model (dot-dashed 
lines) and color glass condensate (dashed lines). The corresponding mass values are m c = 1.5 
GeV and m b = 4.5 GeV. 



complete photon flux instead of the approximation leading to Eq. (12). The final expression 
reads now as 



d ^ Q m = j_m^j fdkURic(k ± ){i + 4 ^-^) arcth * 

dY du 2tt 2 J ± A V ±J \ h,Ju2_L. Am* .F&~ 



where we define the rapidity Y = \n(l/x) = ln(2 uj jl/ Am ,q) and the variable transformation 
between u and Y should be carried out. The equivalent photon flux is taken from Eq. (2). 



4 Results and discussions 



In this section we present the numerical calculation of the rapidity Y distribution and total cross 
section for charm and bottom photonuclear production. In particular, we are focusing mostly 
on LHC domain where small values of x would be probed. At RHIC, x = (Mqq/200 GeV)e~~ y , 
which implies x > 10~ 2 and, consequently, small deviations between the high energy QCD 
approaches. In the following, one considers the charm and bottom masses m c =1.5 GeV and 
m& = 4.5 GeV , respectively. Moreover, for PbPb (A = 208) collisions at LHC, one has the 
c.m.s. energy of the ion-ion system \ASnn = 5500 GeV and the Lorentz factor j L = 2930. 

In Fig. 1 are shown the rapidity distribution, for the distinct high energy approaches con- 
sidered before. The collinear result is denoted by the long-dashed curves, where has been 
used Eqs. (3) and (4) employing the EKS98 parameterization for the collinear nuclear gluon 
function. The solid and dotted lines label the semihard (fcj_-factorization) results, where one 
has used Eq. (6) and the ansatz given by Eq. (7) for the unintegrated gluon function. Two 
possibilities for the nucleon gluon distribution are considered: (I) GRV94(LO) - solid line - and 
(II) GRV98(LO) - dotted line. The saturation model results are denoted by the dot-dashed 
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QQ 


Collinear 


SAT-MOD 


SEMIHARD I (II) 


CGC 


cc 


2056 mb 


862 mb 


2079 (1679.3) mb 


633 mb 


bb 


20.1 mb 


10.75 mb 


18 (15.5) mb 


8.9 mb 



Table 1: The photonuclear heavy quark total cross sections for ultraperipheral heavy ion colli- 
sions at LHC (t/S^ = 5500 GeV) for PbPb. 

lines, with the input given by Eq. (8). The color glass condensate prediction (dashed lines), 
is given by our phenomenological ansatz using expressions (17) and (18). The predictions for 
the collinear approach and the semihard formalism are similar for both charm and bottom 
production and give somewhat larger values than the saturation and CGC results. One possi- 
ble interpretation for the similarity between the predictions of the semihard approach and the 
collinear one is that the expected enhancement in the fc^-factorization formalism, associated 
to the resummation of the (a s In ^|) n in the coefficient function [22], is not sizeable for in- 
clusive quantities in the kinematic region of the future colliders. This feature holds the trend 
already verified at nucleon level as in photon-proton and photon- nuclei collisions [12]. Even 
inclusive cross sections in hadroproduction provide similar results between collinear and semi- 
hard (dipole) approaches [59]. Probably, a more promising quantity to clarify this issue would 
be the transverse momentum p ± distribution. In this case, the semihard approach seems to 
be in better agreement with experimental data in the pp collisions than the collinear approach 
[34]. 

Our phenomenological ansatz within the CGC formalism gives similar results as the satura- 
tion model, but should be noticed that the physical assumptions in those models are distinct. 
While Eq. (9) considers multiple scattering on single nucleons, our expression for the dipole- 
nucleus cross section [Eq.(16)] assumes scattering on a black area filled by partons coming from 
many nucleons. It is important to emphasize that the current experimental data for the nu- 
clear structure function can only be described if the first choice is implemented [39]. However, 
the correct expression for a^ ip in the kinematical range of the future colliders is still an open 
question. 

Let us now compute the integrated cross section considering the distinct models. The results 
are presented in Table 1 for charm and bottom pair production. The collinear approach gives a 
larger rate, followed by the semihard approach, a clear trend from the distribution on rapidity. 
The saturation model and CGC formalisms give similar results, including a closer ratio for 
charm to bottom production. Concerning the CGC approach, our phenomenological educated 
guess for the color field correlator seems to produce quite reliable estimates. 

Lets estimate the uncertainties present in our predictions using the collinear approach. We 
have that for bottom production our prediction decrease by ~ 20 % if we assume mj = 4.75 GeV 
and by ps 10 % if we assume that factorization scale is fi 2 = m 2 .. Moreover, if the MRST gluon 
density [60] is used instead of GRV98 parameterization, our predictions decrease by ~ 10 %. For 
charm production, the differences are larger due to the small values of x probed in that process. 
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Our results for the total charm production cross section in the collinear approach are similar 
to those ones computed in Ref. [21]. However, they differ largely for bottom production, even 
using the same set of scales and parton distributions. We believe that our results are reliable, 
since our prediction for photon-nucleon interaction is consistent with the HERA data [12], as 
well as with the simple expectation cr^A^blx ~ ^ x a -yp->bbx- 

Regarding the semihard approach we have checked the uncertainties coming from quark mass 
and different choices for the gluon parameterization as input for the unintegrated function. In 
comparison with the default value m c = 1.5 GeV for charm, we have an enhancement of ~ 35 % 
using m c = 1.2 GeV, whereas the result decreases the same amount for m c = 1.8 GeV (similar 
results stand for bottom). Considering the default value for the quark mass, the uncertainty 
when using the GRV98 parameterization is of order 20 %. In Table 1 we provide the cross 
section using (I) GRV94 and (II) GRV98 parameterizations. There is an additional uncertainty 
coming from the energy scale /i 2 entering on the strong coupling constant. Here we have used 
the optimal choice, giving correct results at the nucleon level and allowing simple translation 
to the dipole (position space) representation, since the energy scale does not depend on the 
quark transverse momentum. 

In order to check if the differences between the saturation approach and CGC come from 
the integration weights in Eqs. (6) and (18) or from a different mass number dependence for 
the saturation scales, we have used our ansatz Eq. (16) in the calculation of the photonuclear 
cross section. We have found that for lead the CGC ansatz gives a result ~ 40 % lower than 
the saturation model, whereas for calcium the results are almost identical. Therefore, as the 
deviation between saturation model and CGC comes mostly from the different definition for 
the A-dependence of Q s a(x), an experimental analysis of this process for different nuclei can 
be useful. 

In conclusion, the photonuclear production of heavy quarks allow us to constraint already in 
the current nuclear accelerators the QCD dynamics since the main features from photon-nuclei 
collisions hold in the coherent ultraperipheral reactions. As a summary, we have computed the 
photonuclear production of heavy quarks in ultraperipheral heavy ion collisions. One obtains 
the integrated cross section and the rapidity distribution through well established QCD ap- 
proaches, namely the collinear and semihard factorization formalisms as well as the saturation 
model. For the first time, quantitative predictions for the latter two approaches are presented, 
whereas previous collinear calculations are consistently corroborated. Moreover, the color glass 
condensate formalism has been considered using a simple educated guess for the color field 
correlator in the medium, which allowed us reach to reliable estimates at LHC energies. 
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